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PREFACE 


This  report  presents  results  of  a theoretical  study  which  evaluates 
the  response  of  hardened  underground  facilities  in  rock  to  dynamically 
applied  loads  produced  by  explosions.  The  salient  features  of  the  theo- 
retical model  are  summarized  in  charts  suitable  for  preliminary  design 
of  deep  underground  structures.  This  research  was  conducted  by  personnel 
of  the  Phenomenology  and  Effects  Division  (PED)  of  the  Weapons  Effects 
Laboratory  (WEL),  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES), 
during  the  period  January  1975-June  1976. 

The  primary  analytical  development  was  sponsored  by  the  Defense 
Nuclear  Agency  under  Subtask  J3UCAXSX311,  "Underground  Structures 
Studies,"  under  the  guidance  of  Dr.  Kent  Goering.  Computer  parameter 
studies  and  the  preparation  of  this  report  were  sponsored  by  the  Office, 
Chief  of  Engineers  (OCE),  under  Subtask  UA762719ATLo/A1/017 , "Stability 
of  Deep  Underground  Structures  in  Rock,"  which  was  monitored  by 
Mr.  D.  S.  Reynolds. 

This  report  was  written  by  Mr.  J.  R.  Britt,  PED,  under  the  general 
supervision  of  Messrs.  L.  F.  Ingram,  Chief,  PED,  and  W.  J.  Flathau, 

Chief,  WEL. 


COL  G.  [1.  Hilt,  CE,  and  COL  John  L.  Cannon,  CE,  were  Directors  of 
WES  during  the  study  and  preparation  and  publication  of  this  report. 
Mr.  F.  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units 

of  measurement  used 

in  this  report  can  be  con- 

verted  to  metric  (SI) 

units  as  follows: 

Multiply 

By 

To  Obtain 

inches 

2.5h 

centimetres 

pounds  (force)  per 
square  inch 

6.89^*757 

kilopascals 

kilobars 

100.0 

megapascals 

pounds  (force)  per  ^ 
square  inch-second 

6.891*757 

2 

kilopascals-second 

pounds  (force)  per  ^ 
square  inch-second 
per  square  inch 

1.06869  X lo"^ 

kilograms  per  cubic  metre 

degrees  (angle) 

0.0171*5329 

radians 

CHARTS  FOR  PRELIMINARY  DESIGN  OF  DEEP  UNDERGROUND 


STRUCTURES  SUBJECTED  TO  DYNAMIC  LOADS 


CHAPTER  1 
INTRODUCTION 


1 . 1 BACKGROUND 

Design  of  deep  underground  structures  to  resist  the  shock  levels 
induced  by  nuclear  weapons  must  take  into  account  the  strength  of  the 
surrounding  rock.  Hence,  conventional,  dynamic  structural  design  pro- 
cedures^ are  not  adequate  for  designing  these  facilities,  and  design 

2 

engineers  have  resorted  to  static  methods  such  as  outlined  by  Newmark 

3 ^ 

and  Hendron  and  Aiyer.  Drake  and  Britt  have  extended  these  static  pro- 
cedures to  accommodate  dynamic  loads.  The  present  report  summarizes  the 
salient  features  of  this  dynamic  design  method  in  charts  suitable  for 
preliminary  design  of  deep  underground  structures  in  rock. 

1.2  APPLICATIONS  AND  LIMITATIONS  OF  THE  DESIGN  CHARTS 

The  charts  of  this  report  are  applicable  to  the  preliminary  design 
of  cylindrical,  lined  tunnels  in  rock  subjected  to  the  long-duration 
ground  shock  loadings  produced  by  nuclear  weapons.  (it  is  suggested 
that  the  final  design  be  confirmed  by  more  sophisticated  techniques  such 
as  nonlinear  finite  element  computer  codes.)  The  theoretical  model 


Headquarters,  Department  of  the  Army;  "Engineering  and  Design:  Design 

of  Structures  to  Resist  the  Effects  of  Atomic  Weapons";  Technical 
2 Manual  TM  5-856-3,  15  March  1957;  Washington,  D.  C. 

N.  M.  Newmark;  "Design  of  Rock  Silo  and  Rock  Cavity  Linings";  Technical 
Report  TR  TO-llU,  August  1969;  Space  and  Missiles  Systems  Organization, 
^ Air  Force  Systems  Command,  Norton  Air  Force  Base,  Calif. 

A.  J.  Hendron  and  A.  K.  Aiyer;  "Stresses  and  Strains  Around  a Cylin- 
drical Tunnel  in  an  Elasto-Plastic  Material  with  Dllatancy";  Technical 
Report  10,  September  197?;  U.  S.  Army  Engineer  District,  Omaha,  CE, 
Omaha,  Nebr. 

J.  L.  Drake  and  J.  R.  Britt;  "A  Method  for  Designing  Deep  Underground 
Structures  Subjected  to  Dynamic  Loads";  Technical  Report  N-76-9;  Sep- 
tember 1976;  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 


consists  of  multilayered  concentric  cylinders  of  elastoplastic  materials 
with  a time-dependent,  axially  symmetric  load  representing  the  free- 
field  stress  applied  to  the  exterior  boundary.  Each  element  in  the 
cross  section  is  assumed  incompressible  and  its  yield  governed  by  a Mohr- 
Coulomb  failure  criterion.  A first-order  correction  factor  is  given  to 
account  for  compressibility  of  the  materials. 

Experience  has  shown  that  the  present  model  can  be  a useful  design 
tool  if  the  assumptions  of  the  theory  are  met.  In  particular,  the  fol- 
lowing conditions  should  be  approximated  before  using  the  design  charts: 

1.  The  structure  is  composed  of  concentric  cylinders  of  circular 
cross  section  as  shown  in  Figure  1.1. 

2.  The  yield  of  the  structure  and  surrounding  rock  is  governed  by 
the  Mohr-Coulomb  failure  criterion. 

3.  The  expected  deformation  (especially  compaction)  of  liner  and 
rock  is  small. 

U.  The  ultimate  loading  of  the  structure  is  nearly  axisymmetric . 

5.  The  duration  of  the  loading  is  at  least  5 to  10  times  the  re- 
sponse and  engulfment  times  of  the  structure. 

For  example,  use  of  the  charts  appears  to  be  appropriate  for  design  of 
cylindrical  tunnels  with  concrete  and  steel  liners  sited  in  hard  rock 
such  as  granite.  But  the  charts  will  provide  only  a conservative  lower 
bound  estimate  of  the  response  of  structures  located  in  weak,  highly 
compactible  rocks  such  as  tuff. 


CHAPTER  2 


THEORY 

2.1  PROBLEM  FORMULATION 

Solutions  of  the  theoretical  model  are  cast  in  the  form  normally 
used  in  structural  dynamics: 

2 

M ^ = P(t)  - R(e)  (2.1) 

dt"^ 

where : ^ 

M = effective  mass  (mass/length) 
e = diametrical  strain  (positive  inward) 
t = time 

P(t)  = external  load  (positive  in  compression) 

R(e)  = internal  load-resistance 

Definitions  of  the  effective  mass  and  load-resistance  functions  are 
given  in  the  next  section. 

2.2  APPLICATION  OF  THE  ANALYSIS  TO  OBTAIN  DESIGN  CHARTS 

For  the  purposes  of  design  charts,  the  equation  of  motion  (2.1)  and 
its  solutions  are  put  into  a dimensionless  form.  The  effective  resis- 
tance and  mass  are  approximated  with  bilinear  functions  of  normalized 

strain.  Further,  since  the  time  t at  which  the  maximum  strain 

max 

E occurs  for  most  practical  hardened  tunnel-liner  structures  sub- 
max 

jected  to  nuclear  ground  shock  is  small  compared  with  the  duration  of 
the  loading  pulse  or  free-field  stress,  P(t)  is  adequately  approxi- 
mated by  a step  load  with  rise  time  t^  as  shown  in  Figure  1.1.  This 
load  is  a good  approximation  in  most  cases  if  the  duration  of  the  tri- 
angular pulse  describing  the  ground  shock  stress  is  5 to  10  times  t 

md.x 

The  equation  of  motion  (2.1)  is  put  into  dimensionless  form  using 
the  following  elastic  parameters: 


For  convenience,  symbols  and  unusual  abbreviations  used  in  this  report 
are  listed  and  defined  in  the  Notation  (Appendix  B). 


elastic  limit  or  maximum  elastic  strain 


effective  elastic  stiffness 


maximum  elastic  resistance 


effective  mass  (mass/length)  in  elastic  deformation 


The  dimensionless  variables  are  then 


max 


Equation  2.1  becomes 


max 


are  expressed  as  the  bilinear  functions 


R(F)  = 


for  E < 1 


(2.9) 


and 


1 + B(T  - 1)  for  0 > 1 


_ 1 1 for  e 1 

M(e)  = _ _ (2.10) 

(1  + C(e  - 1)  for  e > 1 


(2.11) 


K = — for  e > e (e  > l) 
p Ae  e 


(2.12) 


C 


fe  m 

M Ae 
e 


for 


e > e (e  > l) 
e 


(2.13) 


Examples  of  Bilinear  fits  to  R(e)  and  M(e)  curves  are  shown  in  Fig- 
ure 2.1.  The  lines  1 are  better  approximations  for  e < 2 percent,  but 
lines  2 would  be  used  if  strains  up  to  or  5 percent  were  of  interest. 

The  charts.  Figures  A.l-A.^8,  summarize  the  maximum  values  of  the 

normalized  strain  e = e /e  , and  dimensionless  time  t = 

max  max  e max 

t /T  , at  which  T occurs  for  the  parameters  A , B , and  C . 
max  e max 

Note  that  A is  a measure  of  the  magnitude  of  the  applied  load  or  free- 
field  stress,  B is  a measure  of  the  stiffness  of  the  rock-liner  struc- 
ture in  plastic  deformation,  and  C is  a measure  of  the  change  in  effec- 
tive mass  or  density  as  the  structure  yields.  The  determination  of  B 
and  C and  the  elastic  normalizing  parameters  from  material  properties 
is  discussed  in  the  next  section. 

2.3  DETERMINING  CHART  PARAMETERS  FROM  MATERIAL  PROPERTIES 

The  geometry  of  the  rock-liner  system  is  shown  in  Figure  1.1.  The 

number  of  concentric  material  layers  of  the  system  is  denoted  by  m . 

Generally,  the  liner  is  represented  by  m-1  layers  and  the  rock  by  the 
til 

single  m layer.  Each  material  denoted  by  the  subscript  i is 
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0 1 2 3 

DIAMETRICAL  STRAIN,  % 


Figure  2.1.  Examples  of  bilinear  fits  to  effective 
resistance  and  mass  functions. 


11 


described  by  the  properties  density  , shear  modulus  , cohesion 

, the  angle  of  internal  friction  , and  Poisson's  ratio  . In 

some  cases  the  unconfined  compressive  strength  may  be  known  instead 

of  q . These  parameters  are  related  by 


q = 


1 - sin  <{i 
u 2 cos  (j) 


(2.11*) 


Also,  Young's  modulus  E is  often  given  rather  than  G . Values  of  G 
are  obtained  from  E and  v using 


G = 


2(1  + v) 


(2.15) 


In  Drake  and  Britt's  model,  each  material  layer  yields  from  its 
inner  boundary  r^  outward.  When  a material  layer  is  partially  plastic 
and  partially  elastic,  the  radius  r . denotes  the  boundary  between  the 
interior  plastic  region  and  the  exterior  elastic  region.  These  bound- 
aries are  denoted  by  dashed  lines  in  Figure  1.1.  A material  layer  is 

The  layer  is  completely  plastic  when 

The  effective  re- 


entirely  elastic  when  r . = r. 

pi  1 

r^^  = which  is  the  outer  boundary  of  the  layer, 

sistance  and  mass  functions  R(e)  and  M(e)  can  be  determined  when  the 


values  of  r 


pi 


are  known. 


Drake  and  Britt  give  a relation  for  r^^  involving  acceleration 
or  inertial  terms.  If  the  design  level  of  maximum  diametrical  strain  is 
not  much  larger  than  about  2 to  5 percent,  the  inertial  terms  in  r_ 
can  be  neglected. 

at  the  interior  r 


pi 


A relatively  simple  equation  for  ^ func- 

tion of  diametrical  strain  e at  the  interior  r = r.  of  the  liner  is 

1 


r . 

-Ei  = 

r . 

1 


n.(Y.  + K.  ,e)  + 2q.  v/n.  + 1 

1 1-1  1-1  1 V 1 


l/(n.+2) 


(2.16) 


Drake  and  Britt,  op.  cit. 
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where 


2 sin  A. 

1 

n.  = : — 

1 1 - sin  Ip. 


the  effective  plastic  resistance  of  layer  i is 


n . r n . 

/r  . \ ^ 2q.  \/n”  + 1 /r  . \ ^ 

f =Y  l-Ei)  + —ijLl _ 

i i-l\r.  / n.  \r.  / 

\ 1 / 1 \ 1 / 


and  the  effective  elastic  stiffness  is 


/r  /r  ' /r 

K.  = K.  + 2G.(^  ) 1 -1^  ] 

1 i-l\  r./  i\r . I \r  . I 

\ 1 / \ 1 / \ pi/ 


where  = 0 and  = 0 . In  the  case  n^  = 0 , (ij)^ 


n.-O  ^ \^i 


^ /-El. 


Whenever  Equation  2.l6  gives  a value  of  r ./r.  greater  than  r.  ,/r. 

pi  1 1+1  1 

the  layer  is  completely  plastic  and  r ,/r.  = r.,.,/r.  must  be  used. 

pi  1 1+1  1 

To  calculate  ’ begin  at  i = 1 , the  inside  boundary  of  the 


liner.  In  this  case  Y. 

1-1 


and  K.  . are  zero.  Then  with  r ,/r., 
1-1  d1  1 


known,  Y,  and  K,  can  be  calculated  and  hence  also  r „/r„  . Con- 

11  p2  2 

tinuing  this  process  out  to  the  rock  (i  = m)  surrounding  the  liner,  all 

the  Y.  and  K.  can  be  determined.  In  the  equations  for  Y.  and  K. 
11  11 

the  outside  rock  boundary  Y should  be  taken  as  <»  in  the  case  of 

ra+1 

nuclear  ground  shock  loading.  The  total  resistance  function  R(e)  is 
given  by 


R(e)  = K e + Y 
m m 


. •«  V ■ 


where  Y is  zero  if  the  liner  and  rock  are  elastic,  and  the  effective 
m 

mass  is  obtained  from 


studies  have  shown  that  an  appropriate  value  of  to  use  in  calcu- 

lating M is  r =«  3r  . Equations  2.21  and  2.22  supply  the  effec- 
“ m+1  m 

tive  resistance  and  mass  functions  needed  in  determining  the  parameters 
B and  C of  the  charts. 

To  obtain  B and  C from  R(e)  and  M(e)  , the  elastic  normaliz- 
ing constant  must  first  be  computed  from  the  relations 


where 


^m+1 


«=  and  r,  /r  ->-0  . 
1 m+1 


M = 
e 


(2.23) 


(2.2U) 


where  r 3r  . If  the  inner  layer  of  the  liner  fails  first , then 

m+1  m 

the  elastic  limit  is 


e 

e 


+ 1 


(2.25) 


Using  Equations  2.11,  2.13,  and  2.l6  through  2.25,  values  of  B and  C 
can  be  calculated. 


lU 


In  practice  one  may  not  use  as  given  in  Equation  2.25  but 

rather  fit  R(e)  and  M(e)  to  the  "best"  bilinear  curves.  In  this  way 


e , and  R are  obtained.  Usually, 
e e 


effective  values  of  K , M 

e e 

plotting  a few  points  of  the  R(e)  and  M(e)  curves  is  adequate  to 


allow  one  to  draw  accurate  visual  fits  to  the  curves.  Static  resistance- 

3 

strain  curves,  R(e)  , can  also  be  obtained  using  the  charts  of  Newmark, 

1+ 

and  Hendron  and  Aiyer. 

2.U  A FIRST-ORDER  COMPRESSIBILITY  CORRECTION 

The  equations  described  above  and  the  charts  of  t and  T 

max  max 

(Figures  A.l  through  A.i+8)  are  for  incompressible  media.  The  first- 
order  correction  factor  given  in  Drake  and  Britt ^ may  be  used  to  approxi- 
mately account  for  the  compressibility  of  the  liner  and  surrounding  rock. 

— c 

The  corrected  normalized  maximum  strain,  e , is  related  to  the  incom- 

max 

pressible  value  e”  through 
max 


E = 2(1  - v)e 
max  max 


(2.26) 


where  v is  an  average  of  the  Poisson's  ratios  of  the  liner  and  rock. 
For  common  rock  and  liner  materials,  v ranges  from  about  0.2  to  0.3 
and  2(l  - v)  = 1.5  is  adequate  for  most  preliminary  design 
calculations . 


Newmark,  op.  cit. 

Hendron  and  Aiyer,  op.  cit. 
Drake  and  Britt,  op.  cit. 
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CHAPTER  3 


EXAMPLE  APPLICATIONS  OF  THE  CHARTS 

To  illustrate  the  use  of  the  design  charts,  considea'  the  following 

steel  and  concrete  liner  sited  in  a relatively  weak  granite.  The  inner 

liner  is  steel  which  has  an  inside  diameter  of  U8  inches^  (radius 

= 2U  inches)  and  is  1.5  inches  thick.  Between  the  steel  and  the  granite 

there  is  11  inches  of  concrete.  The  radii  r.  are  then 

1 

r^  = 2h  inches 
r^  = 25.5  inches 
r^  = 36.5  inches 

For  the  resistance  calculation,  take  the  outside  rock  radius  to  be 

r,  = °°  so  that  1/r,  = 0 . For  the  effective  mass  calculation  use 
4 4 

= 3^3  • 

2 

The  material  properties  are 


Steel: 


^ ^ ■ 2,.  2 

= 7*5  X 10  psi-sec  /in 


12  X 10  psi 


= 27,500  psi 


0 deg 


= 0.27 


A table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  3. 

The  units  of  density  were  chosen  to  be  consistent  with  P in  psi 


and  rj  in  inches.  To  convert  unit  weight  in  Ib/ft^  to  density  in 

psi-sec^/in^  divide  by  667,000.  To  convert  density  in  kg/m"^  to  density 

2 2 —8 
in  psi-sec  /in  multiply  by  9-36  x 10 
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Concrete: 


Granite : 


G 


2 

^^2 

^2 


-I4  2,2 

2.25  X 10  psi-dec  /in 

2.13  X 10^  psi 

2,000  psi 
37  deg 

0.25 

-U  2,2 

2.5  X 10  psi-sec  /in 

I4  X lo*’  psi 

2,900  psi 

30  deg 

0.25 


First  determine  the  effective  elastic  stiffness  K , effective 

e 

mass  , and  elastic  limit  using  Equations  2.23,  2.2li,  and  2.2S. 

Expanding  the  summation  in  Equation  2.23  for  three  material  layers, 
obtain 


where  r^/“>  = 0 . Substituting  the  numerical  values  of  the 

parameters  yields 
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I. 


2 

Similarly,  M = 0.231  psi-sec  is  obtained  from  Equation  2.2k.  Using 
0 

Equation  2. IT,  the  value  = 0 is  obtained,  and  from  Equation  2.25 

e = 0.00115  . 
e 

To  determine  resistance  points  in  the  region  of  plastic  deformation 

(e  > e ),  use  Equations  2.l6  through  2.22.  From  Equation  2.17  calculate 
e 

n^  = 0 , n^  = 3.02  , and  n^  = 2.00  . Next,  using  Equation  2.l6  deter- 
mine the  ratios  of  the  plastic  boundaries  r . to  the  inside  radii  r. 

pi  1 

starting  at  i = 1 and  working  up  to  i = 3 . Calculate  for  e = 0.01  . 

When  i=l,  Y.=Y  =0,  K..=K  =0;  hence, 

1-1  o 1-1  o 


- 


UG^e 


l/n^+2 


+ 1 


1(12  X 10°)0.01 


2(27, 500 )V0  + 1 


1/0+2 


= 2.95 


But  ~ 1>0625  . This  means  that  .the  steel  layer  is  completely 

plastic,  and  hence,  ~ ~ 1-0625  must  be  used  for  e = 0.01 

and  larger  strains. 

Yp  and  can  now  be  evaluated  using  Equations  2.l8,  2.19,  and 

2.20.  Equation  2.20  applies  becaxise  n^  = 0 . Since  Y^  = 0 and 
= 0 , the  equations  reduce  to 

r r 

Y.,  = 2q^  \ln,  + 1 In  = 2q.,  In  — = 3331  psi 

1 1 V 1 r^  1 Tp 


and 


These  values  of  Y^  and  do  not  change  for  increasing  strain  e 

since  the  plastic  boundary  r^^  has  reached  its  maximum  value. 

All  information  needed  to  calculate  from  Equation  2.l6  is 

now  available,  and  the  result  is  r _/r„  = 1.33  and  r „ = 33.9  inches 

p2  2 p2 

Thus,  at  the  strain  e = 0.01  the  concrete  Layer  is  about  half  plastic 
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and  half  elastic.  For  i = 2 , Equations  2.18  and  2.19  become 


"2  = ^1  r 


- 1 = 11,520  psi 


/ \ 2 f:  ,2  , ,2' 

K-  = kJ^I  + 2G  (— 1 - [ — 1 

2 1\^2  / ^Vp2j  Ys) 


I 1 \ 5 

- — 1 = 2.93  10^  psi 

V3/ 


These  values  of  and  can  be  substituted  into  Equation  2.l6  to 

obtain  r -,/r_  = 1.15  and  r _ = ^*2.0  inches.  For  i = 3 , Equa- 
p3  3 p3 

tions  2.l8  and  2.19  become,  with  r,  = «>  and  r /r,  = 0 , 


Y =y/j3\  ,!VSZI/!£3\  =16,900  psi 


3 2\r 


/r  \ /r  \ g 

K = KJ-^1  + 2GJ-^j  = 3.00  X 10  psi 

^ "\’'3  / Xp3 


For  the  present  example  the  number  m of  material  layers  is  3. 
Thus  Equation  2.21  for  the  total  resistance  of  the  liner-rock  system 


becomes 


R(e)  = K^e  + 


Substituting  the  values  of  and  Y^  above  and  e = 0.01  yields 

R(O.Ol)  = Ut,000  psi,  which  is  plotted  on  the  lower  graph  of  Figure  2.1. 


From  Equations  2.22  and  2.20  the  effective  mass  M is  calculated: 


M(0.01)  = 0.33  psi-sec^ 


where  the  value 


^1*  = 3r 


3 


has  been  used  as  suggested  previously. 


Using  the  same  procedure  as  for  e = 0.01  , additional  points  at 
e = 0.005  , 0.02  , and  0.03  were  calculated  and  plotted,  along  with 
the  elastic  values,  in  Figure  2.1.  Smooth  curves  were  then  drawn 
through  these  points  and  the  bilinear  fits  denoted  (l)  and  (2)  were  visu- 
ally produced. 


From  previous  experience  with  similar  structures  it  was  determined 

that  the  liner  in  the  example  could  safely  take  about  2 percent  strain 

without  buckling.  This  figure  was  then  taken  as  the  design  criterion  for 

the  structure.  To  apply  the  dynamic  design  charts,  the  bilinear  fits  (l) 

in  Figure  2.1,  which  are  good  up  to  2 to  2.5  percent,  should  be  used. 

From  these  curves  the  effective  elastic  limit  is  e = O.OOU  , the  effec- 

6 ® 

tive  elastic  stiffness  is  ~ ^ 10  psi,  and  the  effective  elastic 

period  is  T = 2it\/M  /K  = 1.1  msec.  These  values  of  e and  K 
e » e e e e 

should  be  used  instead  of  the  results  obtained  earlier  from  Equa- 
tions 2.23  and  2.25.  The  effective  maximum  elastic  resistance  from  the 
fit  (l)  is  then  R = 30,000  psi,  the  effective  plastic  stiffness  is 
Kp  = AR/Ae  ^2.5  X 10  psi,  and  the  slope  of  the  M curve  in  plastic 
deformation  is  AM/Ae  15  . The  chart  parameters  B and  C are  then 

obtained  from  Equations  2.11  and  2.13: 
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^ ^ 2-^  X 10^ 
T.5  X 10^ 


P = _e  M 
^ M Ae 
e 


O.OOU 

0.213 


15  0.26 


The  B = 0.3  curve  of  Figures  A. IT  and  A. lb  should  be  used  to  obtain 


£ = e /e 

max  max  e 


and  t = t /T 
max  max  e 


For  a peak  free-field  stress  P = U3,500  psi  (3  kilcbars), 

^ max  _ 

A = P /R  = 1+3,500/30,000  = 1.1+5  . For  t = 0 Figure  A. IT  gives 
max  e _ r 

T * T and  Figure  A.lb  gives  t =*=  1 . 1 . Hence,  c = e c 
max  max  _ max  e max 

O.OOU  X 7 = 0.02b  or  2.8  percent  and  t = T t =^1.1  msec  x 1.1 

^ max  e max 

1.2  msec.  Using  the  average  value  of  Poisson's  ratio  v = 0.25  , Equa- 
tion 2.26  gives  a compressibility  correction  factor  of  1.5  and  a cor- 
rected maximum  strain  of  = 1+.2  percent.  Thus  for  a suddenly  ap- 

_ max 

plied  load  (t  = O)  and  P = 1+3,500  psi,  the  dynamic  strain  is  more 
r max 

than  double  the  2 percent  design  level  of  maximum  strain  that  has  been 
assumed  for  the  structure. 

Figures  A. 19  thro\igh  A.2U  would  be  used  to  determine  the  structural 
response  if  P(t)  had  a nonzero  rise  time.  In  Figures  A. 23  and  A.21+ 
for  ~ 5 the  response  is  essentially  quasi-static.  These  charts 
show  that  for  the  present  example  the  maximum  static  strain  is  less  than 
half  the  dynamic  strain  produced  by  a suddenly  applied  load. 


At  the  lower  stress  level  P = 29,000  psi  (2  kilobars), 

max 


A = 0.9T  • With  t =0  the  charts  give  values  of  e 

r max 


3 and  t 


0.8  which  yield  e 


1.2  percent  and  t 


V = 0.25  , the  corrected  maximum  strain  is 

max 

within  the  design  level. 


*0.9  msec.  Using 
1.8  percent,  which  is 


1 


APPENDIX  A 


DESIGN  CHARTS 


The  maximum  dynamic  response  of  cylindrical  underground  structures 
subjected  to  the  long-duration  ground  shock  loading  of  nuclear  explo- 
sions is  summarized  in  Figure?  A.l  through  A.U8.  The  maximum  normalized 

diametrical  strain  T = e /s  and  the  dimensionless  time  of  maxi- 

_ max  max  e 

mum  strain  t = t /T  are  given  in  the  charts  for  the  character- 
max  max  e _ 

istic  parameters  A , B , C , and  t where:  i 

r 

T^  = 2ir^M^/K^  = effective  elastic  period 

A = P /R  = dimensionless  maximum  applied  load 
max  e 

B = K /K  = dimensionless  plastic  stiffness 
p e 

C = ( Sg/Me ) ) = dimensionless  rate  of  change  of  effective  mass 

in  plastic  deformation 

t = t /T  = normalized  rise  time  of  applied  load 

The  material  properties  needed  to  determine  the  chart  parameters  are  il- 
lustrated in  Figure  2.1;  an  example  using  the  charts  is  given  in  Chap-  - 

ter  3.  Charts  are  included  in  this  appendix  for  normalized  rise  times  j 

ranging  from  t^  = 0 for  a suddenly  applied  load  to  t^  = 2 or  5 where  | 

the  structural  response  is  essentially  quasi-static.  On  the  figures  for  * 

t = 2 or  5 some  of  the  t curves  for  the  larger  values  of  B have  I 

r max  ^ 

been  smoothed  to  produce  readable  "average"  curves  which  deviate  from  3 

the  original  ones  by  less  than  about  one  fourth.  ■ 
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XDUl  XOUi 


max  max 


max  max 


Fiffure  A. 13*  £ curves  for  tunnel-liner  systems  with  C - 0.1  and 

® max 


max  max 


r ^ 


0 2 4 6 8 10 

1 3 5 7 9 

A = P.ox/R. 


Figure  A.2U.  t curves  for  tunnel-liner  systems  with  C = 0.3 
max 

and  t = 5 • 
r 


i 


J 

J 


I 


i 


Figure  A. 30.  t curves  for  tunnel-liner  systems  with  C = 0.5 
max 

and  t = 1 . 
r 


mox  max' 


Figure  A. 31. 


e curves  for  tunnel-liner  systems  with  C 
max 


0.5  and 


t 

r 


5 . 


max  max 


Figure  A. 3b.  t curves  for  tunnel-liner  systems  with  C = 0.7 
max 

and  t =0.5. 

r 


J 


1 


max  max 


Figure  A.UU. 


P /R 

max  e 


t curves  for  tunnel-liner  systems  with 
max 

and  t =0.5* 
r 


Fiimre  A.U6.  t curves  for  tunnel-liner  systems  with  C = 1 
^ max 

and  t = 1 . 
r 


Figure  A.  ^8.  t curves  for  tunnel-liner  systeins  with  C = 1 
^ max 

and  t = 5 . 
r 


APPENDIX  B 


C 


A = P /K 
max  e 

B = K /K 
P e 

(e  /M  )(AM/Ae) 
e e 


E 

G. 

1 

i 

K 

e 

K. 

1 

K = AR/Ae 
P 

m 


M 

M 

_ e 

M(D  = M/M 

e 

n. 

1 

P 

max 

p(t) 


r . 

1 

r . 

pi 


R = K G 
e e e 

R(e) 

R(r)  = R(g)/R 

e 

t 

t = t/T 

e 

t 

_ max 

t = t /T 
max  max  e 


t 

r 


NOTATION 

Dimensionless  peak  free-field  stress 
Dimensionless  plastic  stiffness 

Dimensionless  rate  of  change  of  M(e')  in  plastic 
def  ormation 

Young's  modulus 

Shear  modulus 

Subscript  referring  to  a component  material  layer 
Effective  elastic  stiffness 
Effective  elastic  stiffness  of  layer  i 
Effective  plastic  stiffness 

Number  of  material  layers  in  liner-rock  system 
Effective  mass  (mass/length) 

Effective  mass  (mass/length)  in  elastic  deformation 
Dimensionless  effective  mass  (mass/length) 

2 sin  4>^/(l  - sin  4>^) 

Peak  value  of  P(t) 

External  load  (or  free-field  stress) 

Cohesion 

Radius  to  inside  of  material  layer  i 

Radius  to  the  outer  boundary  of  the  plastic 
region  of  material  layer  i 

Maximum  elastic  resistance 

Internal  load-resistance 

Dimensionless  resistance 

Time 

Dimensionless  time 

Time  at  which  e occurs 
max 

Dimensionless  time  at  maximum  strain 
Rise  time  of  external  load 
Normalized  rise  time  of  P(t) 

Effect 've  elastic  period 
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Effective  plastic  resistance  of  layer  i 

€ Diametrical  strain,  positive  inward 

T = e/f  Normalized  diametrical  strain 
e 

e Elastic  limit,  maximum  elastic  strain 
e 

e Maximum  value  of  diametrical  strain  e(t) 

_ max 

e = e /e  Normalized  maximum  strain 

max  max  e 

Q 

€ t corrected  for  compressibility 

max  max 



e e corrected  for  compressibility 

max  max 

V.  Poisson's  ratio 

1 

Density  of  material  layer  numbered  i 
Unconfined  compressive  strength 
(Ji.  Angle  of  internal  friction 
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